Temperature dependence of the conductivity of graphene on boron nitride 
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The substrate material of monolayer graphene influences the charge carrier mobility by various 
mechanisms. At room temperature, the scattering of conduction electrons by phonon modes local- 
ized at the substrate surface can severely limit the charge carrier mobility. We here show that for 
substrates made of the piezoelectric hexagonal boron nitride (hBN), in comparison to the widely 
used Si02, this mechanism of remote phonon scattering is -at room temperature- weaker by almost 
an order of magnitude, and causes a resistivity of approximately 3 Q. This makes hBN an excellent 
candidate material for future graphene based electronic devices operating at room temperature. 

PACS numbers: 72.80.Vp,77.84.Bw,72.10.Di 



I. INTRODUCTION 

Apart from the great theoretical interest in the elec- 
tronic properties of graphen^ spurred by the experimen- 
tal realization of these one atom thick carbon layerSj^^l 
the new material has been soon recognized as a promis- 
ing candidate for technological applications.!^ At low tem- 
peratures, electrons in free standing (that is, suspended 
above the substrate) graphene sheets can travel in an 
essentially ballistic manner, without scattering over dis- 
tances of the order of a micron.!^ With the carrier den- 
sity tuned to finite values by an applied gate voltage, the 
temperature dependence of the resistivity in suspended 
graphene sheets shows a metallic behavior: at temper- 
atures above 50 K, it increases linearly with tempera- 
ture. The carrier mobility at room temperature in sus- 
pended graphene is severeh' limited by scattering by out 
of plane flexural phonons,"^ resulting in a resistivity that 
increases quadratically with the temperature. The effect 
of flexural phonons on the charge carrier mobility can be 
eliminated by applying strain or placing graphene on a 
substrate. 13 

Compared to graphene devices on substrates, sus- 
pended graphene samples offer charge carrier mobilities 
which are higher by a factor of ten, approximately."^ Like- 
wise, the thermal conductivity of graphene on a substrate 
can be an order of magnitude lower than that of sus- 
pended graphene, due to damping of the flexural phonon 
modes.lS'ii' For technological applications, however, the 
suspended geometry imposes rather strict limitations on 
the device architecture. Moreover, free floating graphene 
is always crumpled, with the corrugations inducing ef- 
fective fields which in turn influence the charge carriers. ^'^ 
While many experiments so far involved graphene sam- 
ples deposited over Si02 or grown over SiC substrates, 
hexagonal boron nitride (hBN) has developed into a 
promising candidate as a sub strate material for improved 
graphene based devices. 1^^^^^ hBN is a piezoelectric, large 
band gap insulator isomorphic to graphite. Boron and ni- 
trogen atoms occupy the inequivalent sublattices in the 
Bernal structure. The lattice mismatch with graphite is 



small (1.7%),^^ and hBN is expected to be free of dan- 
gling bonds and surface charge traps. 



A number of mechanisms limit the low temperature 
carrier mobility in graphene .l^^l^ At room temperature, 
interaction with optical surface phonon modes on the 
interface between graphene and the Si02 substrate was 
found to play an important role.l^Sl It is this mechanism 
of remote-phonon scattering from an hBN substrate that 
we want to consider in the present work. In general, the 
polar phonon modes on the surface of the substrate (with 
energies of 50—200 meV) create a long-range electric field, 
which influences the electrons in the graphene sheet, typ- 
ically around 4 A away. The influence of remote phonon 
scattering on the carrier mobility in two-dimensional 
electron systems is a well-known phenomenon in semi- 
conductor physics, and was investigated for quantum 
wells and other heterostructures including metal-oxide- 
semiconductor field-effect transistors (MOSFETs).'2iJ22| 
The effect is more pronounced in graphene due to the 
much smaller vertical dimension of the devices, as deter- 
mined by the van der Waals distance. Also, the band gap 
in semiconductor systems prevents low energy interband 
transitions, which are present in graphene with the two 
bands touching at the Dirac point or with the metallic 
band in the case of doping. For single-layer graphene 
sheets on substrates of SiC or Si02, surface-phonon scat- 
tering has been investigated in Refs. employing 
various methods for calculating the conductivity and de- 
scribing the screening of the interaction by conduction 
electrons. In the present paper we calculate the resistiv- 
ity due to surface phonon scattering for hBN substrates. 
We show below that the temperature dependence of this 
mechanism scales with the thermal population of surface 
phonon modes. While for Si02 substrates, the resistivity 
due to remote phonon scattering is known to be com- 
parable or might surpass that due to graphene intrinsic 
phonons, we show that for hBN, the effect is almost 
an order of magnitude smaller, resulting in desirable high 
charge carrier mobilities. 
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II. ELECTRON SELF ENERGY DUE TO 
REMOTE PHONON COUPLING 

A. Surface phonons 

The dielectric function of a substrate material with 
transverse optical modes reads 
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where eoo denotes the high frequency dielectric constant 
of the material and the diniensionless oscillator strengths 
fi measure the contribution of each mode to the screen- 
ing properties of the material. They are determined 
from experimental data by defining intermediate dielec- 
tric constants e^. These are evaluated at a frequency just 
above the corresponding resonance (see for exam- 
ple Ref. m]). With eo the static dielectric constant, the 
oscillator strengths are given by fi = e^-i — e^. The 
frequencies Wg of the corres pondin g surface modes are 
determined from the equatiorPHSl! 

€{0J) + 1=0, 

where 1 is the dielectric constant of air and we neglected 
the dielectric response of the (atomically thin) graphene 
layer. 

The polar surface modes on the substrate interface cre- 
ate a polarization field which decays exponentially with 
the distance from the interface, and is felt by the elec- 
trons in the graphene sheet. This remote interaction can 
be brought to the form^^ 



i k,q 



where i runs over the different surface modes ujq, k and 
q are two-dimensional momentum vectors parallel to the 
graphene-substrate interface, Cq and c!!q are destruction 
and creation operators for surface phonons and the 
and fek operators are the destruction operators for elec- 
trons on the A and B sublattices of the graphene sheet. 
They are coupled by the interaction matrix element 



(1) 



where z denotes the (positive) distance between the sub- 
strate and the graphene sheet, the coupling strength to 
the individual surface modes is given by 
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+ 1 e»-i + 1 
and the dimensionless function reads 
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TABLE L Material parameters for hBN, Si02 and SiC sub- 
strates, taken from Refs. ''2F, 1211 and 1291 respectively. The 
bulk optical phonon frequencies ljto are given at the F point. 
For materials with two surface modes, eoo figures as the inter- 
mediate dielectric constant £2 in the calculation of the cou- 
pling parameter 7^ in Eqn. |2|. 
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eo 


5.09 


3.90 


9.7 


ei 


4.57 


3.05 




^00 


4.10 


2.50 


6.5 


hjj^o [meV] 


97.4 


55.6 


97.1 


huj^o [meV] 


187.9 


138.1 




7^ 


0.0153 


0.0428 


0.040 


7^ 


0.0165 


0.0388 




hioo [meV] 


101.6 


61.0 


116 


hujo [meV] 


195.7 


149.0 




z[A] 


3.4 


4.0 


4.0 



with A denoting the surface area of the interface, e 
the electron charge, and Cvac the permittivity of free 
space. The inverse Thomas-Fermi screening length 
qtf = e^^'F/(7revac^^^'|-) incorporates the effect of dy- 
namic screening from the conduction electrons. 



B. Electron self energy 

In considering the charge carriers in graphene, we limit 
our treatment to the Dirac cone approximation, where 
the dispersion reads 



svf\P\ 



(4) 



with s = 1 denoting the conduction (tt*) band and s — 
— 1 the valence (tt) band, and the Fermi velocity vp ^ 
10^ m/s. In leading order perturbation theory, the self 
energy acquired by these Dirac fermions due to coupling 
to a remote substrate phonon is 

S,(p,ip)= p ^ s 1^ p^^.p (5) 
= -fcsT J2 M'i?^°'(«^n)Gi")(p + q,ip-t-fw„) , 

where D^^^ and G^^-' denote the free thermal (Matsub- 
ara) Green's function for the phonon and the electron, 
respectively, and the interaction ([T]) is depicted by the 
vertical zigzag lines in the diagram. (See Eqns. ( Al ) and 
(A3) in the Appendix and Refs. [30l and [311 for details.) 
After summation over the bosonic Matsubara frequency 
LOn , a shift of the integration variable to k = p + q, and 
rotating back to real frequencies {ip uj — ^ + iif), we 
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are left with 



h{uj — /i — lo'q) — Eg{k) + irj 



N'^ + nF[E,{k)] 



h{uj — /i + Wq) — Es{k) + irj 



g,{\.){Mi_^f (6) 



Nl the ther- 



Here, np denotes the Fermi distribution, 
mal occupation of the phonon, and the appearance of 
the 2x2 matrix Qg in the electron Green's function (see 
Eqn. ( A4 )) is due to the spinor representation of the elec- 
tron wavef unction in graphene. 



III. QUASIPARTICLE SCATTERING RATE 

The scattering rate can be obtained from Eqn. (|6| by 
multiplying the self energy from the left and right with 
spinor wavefunctions F±i (see Eqn. ( A2 )) for the ingoing 
and outcomming electron, respectively, and afterwards 
taking the imaginary part: 



2h 



=±1 
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where the incoming momentum p is set on-shell, that is 
h(jj — zLhvF\p\ — and sgn[a;] = ±1 for w > (w < 0), 
respectively. The imaginary part is obtained via lm[x + 
ir]]~^ = —ttS{x), which yields the following rates T±{(jj) 
for emission or absorption of a phonon with frequency 

, ,i . 



(8a) 
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for < w < Wq 
else 



T\iu;) = ^[No + l-nF{io~iol,)] 
Y,M^_^SivFk-\u:-ul\)l^^-^' 

^-i^) = ^[No + nFiio + ,,i,)] (8b) 

E, .9 f / , I , IN I /-I for — < w < 

M^-kSivFk-\oj+oj^\)r^^' g^gg °- 

The total scattering rate is obtained by summing the ab- 
sorption and emission rates for scattering from all surface 
phonons i. The angular factor 

/±i(k,p) = ^[l±cos(0k-0p)] , 

in Eqns. where k = fc(cos6'k, sin^k), distinguishes 
between intraband (upper sign ) and interband (lower 
sign) scatterin^2ll^ gge Eqn. ( A5 ) and the inset in Fig. [l] 
Terms describing interband scattering are seen to con- 
tribute only in the range |a;| < ujq. 

Figure fl] shows r(a;) at T = 300 K for intrinsic 
graphene [Ep = 0, full line) and extrinsic graphene 
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FIG. 1. Quasiparticle scattering rate r(cj) in single layer 
graphene due to surface phonons of the hBN substrate at 
T — 300 K [see eqns. Q]. The frequencies of the optical sur- 
face modes uig''^ are given in table |lj the distance between 
graphene layer and substrate is 3.40 A. Full line: Intrinsic 
graphene (Ep = 0). Dashed line: Extrinsic graphene with 
Ep = 0.5 eV. Inset: The angular factor f±i [see Eqn. (A5l] 



in the scattering rate distinguishes between interband and in- 
traband scattering. 



with Ep — 0.5 eV (corresponding to a carrier density 
of 71 = 1.83 • 10^'^ cm^^, dashed line) on a hBN substrate. 
The distance between graphene layer and substrate is set 
to 3.40 A, as found in Ref . fT6l for a stacking configuration 
with one carbon over N, and the other carbon centered 
above a hBN hexagorP^l. For electron energies around 
Ep ± hujQ , scattering is strongly suppressed and goes to 
zero at T = 0, as there are no empty electronic states 
below the Fermi energy to scatter into. 

The large difference in F between the doped and un- 
doped case is due to the Thomas-Fermi screening in 
Eqn. (|3|. As this model assumes the instantaneous re- 
action of the screening charges, our values for F at fi- 
nite doping present a lower bound on the scattering 
rate.'23Elljf screening is completely neglected, the rate at 
Ep = 0.1 eV is larger by a factor of four, approximately, 
which presents an upper bound for the rate. 



IV. SUBSTRATE LIMITED CONDUCTIVITY 

We calculate the electrical dc conductivity a of the 
graphene layer via the Boltzmann equation: 



e 



cLu 



Ftr(t^) 



dnp{uj) 
dco 



(9) 



Here, the transport scattering rate Ftr is defined as in 
Eqn. ([8|, but with an additional factor [1 ^ cos(6'k — dp)] 
in the integrand (with the upper sign for intraband and 
the lower one for interband scattering) , which lends more 
weight to large angle scattering events. 1^ As the deriva- 
tive of the Fermi function is sharply peaked around the 
Fermi energy, the integrand of Eqn. ^ depends mainly 
on the scattering rate for ui « Ep/h, where both the 
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FIG. 2. Scattering rate r{uj) at T = 300K, Ep = 0.1 eV, 
for different substrate materials. The plot shows the energy 
range around the Fermi energy, which is relevant for the con- 
ductivity, for hBn, SiC and Si02 substrates. (See table |l] for 
the material parameters.) The full line is evaluated with the 
distance between hBN and graphene set to 3.4 A (see Ref. ll6p . 
For comparison, the dots show the rate at a distance of 4 A, 
as for the other two materials (see Ref. I24p . The resulting 
relative difference in V{Ef) is less than 5 per cent. 
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FIG. 3. Temperature dependence of the resistivity for Si02, 
hBN and SiC (top to bottom). The dots corresp ond to the 
population factor p()/ s'mh[hu)li / (ksT)] (see Eqn. (lOl), with 
po = 178,95,381^ for Si02, hBN, and SiC, respectively. 
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FIG. 4. Conductivity versus Fermi energy for different tem- 
peratures T = 300, 275, 250 K (bottom to top). 



phonon emission and phonon absorption term in Eqn. ([8]) 
are proportional to 



ul/Bin\,[hul/{kBT)] 



(10) 



At a given temperature, this is a decaying function of 
Wq- Thus scattering in the relevant range is larger 
for substrate materials with lower phonon frequencies, 
which translates into a lower conductivity. In Fig. [2] 
we compare scattering rates at electron energies near 
Ep — 0.1 eV for different substrate materials. For elec- 
trons at the Fermi energy, the rates for SiC and hBN 
differ roughly by a factor of two, while the scattering 
rate on Si02 is around ten times larger than that on SiC. 
As the resistivity is calculated from Eqn. ([9]), its temper- 
ature dependence, shown in Fig. |3j is likewise dominated 
by the population factor (10). We here assume a tem- 



perature which is constant throughout the sample, with 
local heating of the graphene layer prevented by heat dis- 
sipation through the substrate.'^S! 

Regarding the dependence on the charge carrier con- 
centration, the conductivity ^ depends approximately 
linear on Ep for energies larger than Wq (see Fig. |4|. 
This behaviour is intermediate between that typical for 
short range scattering, a{Ep) « const., and for charged 
impurity scattering, a{EF) ss EpW^ 



V. SUMMARY AND CONCLUSIONS 

Wc have analyzed the scattering of graphene electrons 
from phonons localized at the substrate surface. The 
temperature dependence of the scattering rate for elec- 
trons with energies near Ep is proportional to the ther- 
mal population of these surface phonons, and the same 
holds for the resistivity. The resistivity induced by the 
polar surface phonons of hBN substrates at room tem- 
perature is of the order of 3f2, an order of magnitude 
smaller than the corresponding resistivity of graphene on 
Si02 substrates, in consistency with the higher frequency 
of these modes in hBN. 

For graphene - Si02 devices, the temperature de- 
pendence of the conductivity was experimentally found 
to consist of two contributions P2l scattering by acoustic 
phonons in graphene and scattering by surface phonons of 
the substrate. The former shows a linear T dependence, 
is independent of the carrier density, and contributes ap- 
proximately 30 at room temperature. The latter is 
carrier density dependent, follows the bosonic popula- 
tion of the surface modes, and becomes only relevant for 
T > 200 K. At room temperature it is found to dominate 
the linear term. In experiments with monolayer graphene 
on hBNjtS! the conductivity shows the same linear tem- 
perature dependence as that reported in Ref.UDlfor Si02. 
No indication of activated remote surface phonon scat- 
tering was seen up to the experimentally realized temper- 
atures of 200 K. An estimate of the temperature depen- 
dence of the resistivity of graphene on hBN can also be 
obtained from the experiments reported in Ref. 15, The 
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result is consistent with the assumption that the room 
temperature resistivity of graphene on hBN is mainly de- 
termined by in plane phonons. 

Apart from electron-phonon scattering, Coulomb scat- 
tering by charged impurities gives rise to a temperature 
independent residual resistivity. At carrier concentra- 
tions of lO^^cm"^, this is reported with approximately 
400 H. for S i02 and roughly three times smaller values for 
hBN™ 



Our calculation suggests that for graphene on hBN, 
the temperature-dependent part of the resistivity is -in 
contrast to Si02 substrates- at room temperature still 
dominated by electron scattering from the graphene in- 
trinsic phonons, the contribution from interface phonons 
being negligible. Together with the smaller residual re- 
sistivity of graphene on hBN as compared to Si02, this 
might allow for graphene - hBN devices with charge car- 
rier mobilities close to that of suspended graphene. 



Appendix A: Matsubara Green's functions and 
electronic wavefunctions 



The Fourier transform of the free thermal Green's func- 
tion for a surface phonon with frequency ujq is given by 



2hujn 



(Al) 



where the bosonic Matsubara frequencies are defined as 
ujn = 2ti k bT n with integer 

The wavefunction for electronic states in graphene near 
one of the Dirac points is^ '^^ 

F,(r) = A-i/2F^(q)exp(iq.r) 

where s = ±1 denotes the band index, A is the area of 
the system and 



F.(q) = 



e 

2 \ s 



(A2) 



In the same spinor representation, the electron Green's 
function is written as the 2x2 matrixes 

G(°)(k,*c.„) = g,{\^), \ (A3) 

where a;„ — {2n + \)'KkBT are the fermionic Matsubara 
frequencies. Eg is the energy of the electron (see Eqn. Q) 
within the Dirac cone approximation and 
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The angular factor f±i appearing in Eqn. ([8]) is composed 
of 



1, 



FT,(p)g,(k)F,,(p) = -[l + ss'cos(0k-^p)] = /s.s'(k,p) 

(A5) 
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